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Abstract 

Weierstrass's everywhere continuous but nowhere differentiable function is shown to be 
locally continuously fractionally differentiable everywhere for all orders below the 'critical 
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^ | ■ order' 2 — s and not so for orders between 2 — s and 1, where s, 1 < s < 2 is the box dimension 

] of the graph of the function. This observation is consolidated in the general result showing 

a direct connection between local fractional differentiability and the box dimension/ local 
Holder exponent. Levy index for one dimensional Levy flights is shown to be the critical order 
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of its characteristic function. Local fractional derivatives of multifractal signals (non-random 
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. derivatives provide a powerful tool to analyze pointwise behavior of irregular signals. 
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functions) are shown to provide the local Holder exponent. It is argued that Local fractional 
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Attractors of some dynamical systems are examples of the occurrence of continuous 
but highly irregular (nondifferentiable) curves and surfaces. Frequently their graphs 
are fractal sets. Ordinary calculus is inadequate to characterize and handle such curves 
and surfaces. In this paper we evolve the notion of "local fractional derivative" by 
suitably modifying the concepts from fractional calculus, a branch which allows one 
to deal with derivatives and integrals of fractional order. In particular we establish a 
direct quantitative connection between the local scaling behaviour (or dimension) and 
the order of differentiability. The bigger the fractal dimension the smaller is the extent 
of differentiability. We show that the method developed here provides a powerful tool 
for analysis of irregular and chaotic signals. It is further noted to be suitable for dealing 
with fractal processes. We have also established a local fractional Taylor expansion, 
which should be of value in the approximation of scaling signals and functions. 



1 Introduction 



The importance of studying continuous but nowhere differentiable functions was emphasized a long 
time ago by Perrin, Poincare and others (see Refs. and 0). It is possible for a continuous 
function to be sufficiently irregular so that its graph is a fractal. This observation points out 
to a connection between the lack of differentiability of such a function and the dimension of its 
graph. Quantitatively one would like to convert the question concerning the lack of differentiability 
into one concerning the amount of loss of differentiability. In other words, one would look at 
derivatives of fractional order rather than only those of integral order and relate them to dimensions. 
Indeed some recent papers 0, f|, |[ || indicate a connection between fractional calculus |7|, [| and 
fractal structure 



T0|| or fractal processes [XT], [L2], |13| . Mandelbrot and Van Ness [[TI 



have used 
In Refs. 
Also 



fractional integrals to formulate fractal processes such as fractional Brownian motion. 
J|] and || a fractional diffusion equation has been proposed for the diffusion on fractals. 
Glockle and Nonnenmacher have formulated fractional differential equations for some relaxation 
processes which are essentially fractal time |U| processes. Recently Zaslavasky (T^| showed that 
the Hamiltonian chaotic dynamics of particles can be described by a fractional generalization of 
the Fokker-Plank-Kolmogorov equation which is defined by two fractional critical exponents 
responsible for the space and time derivatives of the distribution function correspondingly. However, 
to our knowledge, the precise nature of the connection between the dimension of the graph of a 
fractal curve and fractional differentiability properties has not been established. 

Irregular functions arise naturally in various branches of physics. It is well known that the 
graphs of projections of Brownian paths are nowhere differentiable and have dimension 3/2. A 
generalization of Brownian motion called fractional Brownian motion |2|, [K| gives rise to graphs 
having dimension between 1 and 2. Typical Feynmann paths [[15, TI|, like the Brownian paths are 



continuous but nowhere differentiable. Also, passive scalars advected by a turbulent fluid |T7|, |TS[ 
can have isoscalar surfaces which are highly irregular, in the limit of the diffusion constant going 



to zero. Attractors of some dynamical systems have been shown |19| to be continuous but nowhere 
differentiable. 

All these irregular functions are characterized at every point by a local Holder exponent typically 
lying between and 1. In the case of functions having the same Holder exponent h at every point it 
is well known that the box dimension of its graph is 2 — h. Not all functions have the same exponent 
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h at every point but have a range of Holder exponents. A set {x\h(x) = h} may be a fractal set. 
In such situations the corresponding functions are mult if r act al. These kind of functions also arise 
in various physical situations, for instance, velocity field of a turbulent fluid |20] at low viscosity. 



Also there exists a class of problems where one has to solve a partial differentiable equation 
subject to fractal boundary conditions, e.g. the Laplace equation near a fractal conducting surface. 
As noted in reference [[H] irregular boundaries may appear, down to a certain spatial resolution, to 



be non-differentiable everywhere and/or may exhibit convolutions over many length scales. Keeping 
in view such problems there is a need to characterize pointwise behavior using something which can 
be readily used. 

We consider the Weierstrass function as a prototype example of a function which is continuous 
everywhere but differentiable nowhere and has an exponent which is constant everywhere. One 
form of the Weierstrass function is 

oo 

W x {t) = ]T A (s ~ 2)/c sin A fc t, t real. (1) 

k=l 



For this form, when A > 1 it is well known [22] that W\{t) is nowhere differentiable if 1 < s < 2. This 



curve has been extensively studied [jl], [H| |24| and its graph is known to have a box dimension s, 



for sufficiently large A. Incidently, the Weierstrass functions are not just mathematical curiosities 
but occur at several places. For instance, the graph of this function is known [|II], [19! ^° ^ e a 



repeller or attractor of some dynamical systems. This kind of function can also be recognized as 



the characteristic function of a Levy flight on a one dimensional lattice |25[ , which means that such 
a Levy flight can be considered as a superposition of Weierstrass type functions. This function has 
also been used |Hj to generate a fractional Brownian signal by multiplying every term by a random 
amplitude and randomizing phases of every term. 

The main aim of the present paper is to explore the precise nature of the connection between 
fractional differentiability properties of irregular (non-differentiable) curves and dimensions/ Holder 
exponents of their graphs. A second aim is to provide a possible tool to study pointwise behav- 
ior. The organization of the paper is as follows. In section II we motivate and define what we 
call local fractional differentiability, formally and use a local fractional derivative to formulate the 
Taylor series. Then in section III we apply this definition to a specific example, viz., Weierstrass' 
nowhere differentiable function and show that this function, at every point, is locally fractionally 
differentiable for all orders below 2 — s and it is not so for orders between 2 — s and 1, where s, 
1 < s < 2 is the box dimension of the graph of the function. In section IV we prove a general result 
showing the relation between local fractional differentiability of nowhere differentiable functions 
and the local Holder exponent/ the dimension of its graph. In section V we demonstrate the use 
of the local fractional derivatives (LFD) in unmasking isolated singularities and in the study of the 
pointwise behavior of multifractal functions. In section VI we conclude after pointing out a few 
possible consequences of our results. 
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2 Fractional Differentiability 



We begin by recalling the Riemann-Liouville definition of the fractional integral of a real function, 
which is given by |7|. |j| 

<pf{ x ) i r f(y) 



[d(x-a)]i T(-q)Ja (x-yY +1 
and of the fractional derivative 

d«f(x) 1 d f* f(y) 



J a ( x -l) q +i dy for q < °' a rea1 ' ( 2 ) 



-dy for < g < 1. (3) 



[d{x-a)]i T(l-q)dxJa (x - y)i 
The case of q > 1 is of no relevance in this paper. For future reference we note |7|. ^ 

dw r(p + i) __. , ... 

— — = — — '—x v q for p > -1. (4) 

dxi T{p-q+l) W 

We also note that the fractional derivative has the property (see Ref. [0]), viz., 

d q f(Px) 

which makes it suitable for the study of scaling. 

One may note that except in the case of positive integral q, the qth derivative will be nonlocal 
through its dependence on the lower limit "a". On the other hand we wish to study local scaling 
properties and hence we need to introduce the notion of local fractional differentiability. Secondly 
from Eq. @) it is clear that the fractional derivative of a constant function is not zero. These two 
features play an important role in defining local fractional differentiability. We note that changing 
the lower limit or adding a constant to a function alters the value of the fractional derivative. This 
forces one to choose the lower limit as well as the additive constant before hand. The most natural 
choices are as follows. (1) We subtract, from the function, the value of the function at the point 
where fractional differentiability is to be checked. This makes the value of the function zero at that 
point, washing out the effect of any constant term. (2) The natural choice of a lower limit will be 
that point, where we intend to examine the fractional differentiability, itself. This has an advantage 
in that it preserves local nature of the differentiability property. With these motivations we now 
introduce the following. 

Definition 1 If, for a function f : [0, 1] — > M, the limit 

J ( J *->v d(x -y) q 

exists and is finite, then we say that the local fractional derivative (LFD) of order q, at x = y, 
exists. 

Definition 2 We define critical order a, at y, as 

ct(y) = Sup{q\a\\ local fractional derivatives of order less than q exist at y}. 
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Incidentally we note that Hilfer PU| , [Z7J] used a similar notion to extend Ehrenfest's classification of 
phase transition to continuous transitions. However in his work only the singular part of the free 
energy was considered. So the first of the above mentioned condition was automatically satisfied. 
Also no lower limit of fractional derivative was considered and by default it was taken as zero. 

In order to see the information contained in the LFD we consider the fractional Taylor's series 
with a remainder term for a real function /. Let 

f ,„_ n)= *Wi (7) 
(V ' V ' q> [d(x-y)]- ■ ( > 



It is clear that 



Now, for < q < 1, 



E>«f(y) = F(y,0;q). (8) 



t(x) - /(») - ^rjg*^ <m 

- J-[F(y,t;q) I \x - y -t^dt]^ 



4- 

T(q) Jo dt 



1 pdF(y,t; g )(x-y-t)^ (1Q) 



provided the last term exists. Thus 



f(*)-f(y) = ^T^-y) 11 



Jo -^r—(*-y-t) q dt, (ii) 



r( 9 + i) Jo 

i.e. 



/(*) = m + Y(ivi) {x ~ y)q + Ri{x > yl (12) 

where Ri(x,y) is a remainder given by 

^ 1 [*-v dF(y,t;q) f 

Rl ^ y) = T{q-TT)Jo -^^(x-y-t) 9 dt (13) 

Equation ( |T3 ) is a fractional Taylor expansion of /(x) involving only the lowest and the second 
leading terms. This expansion can be carried to higher orders provided the corresponding remainder 
term is well defined. 

We note that the local fractional derivative as defined above (not just the fractional derivative) 
provides the coefficient A in the approximation of f(x) by the function f(y) + A(x — y) q /T(q+ 1), for 
< q < 1, in the vicinity of y. We further note that the terms on the RHS of Eq. (O) are non-trivial 
and finite only in the case q = a. Osier in Ref.[^] has constructed a fractional Taylor series using 



usual (not local in the present sense) fractional derivatives. His results are, however, applicable 
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to analytic functions and cannot be used for non-differentiable scaling functions directly. Further 
Osier's formulation involves terms with negative q also and hence is not suitable for approximating 
schemes. 

One may further notice that when q is set equal to one in the above approximation one gets the 
equation of the tangent. It may be recalled that all the curves passing through a point y and having 
the same tangent form an equivalence class (which is modeled by a linear behavior). Analogously 
all the functions (curves) with the same critical order a and the same D a will form an equivalence 
class modeled by x a [ If / differs from x a by a logarithmic correction then terms on RHS of Eq. 
(|TT|) do not make sense precisely as in the case of ordinary calculus] . This is how one may generalize 
the geometric interpretation of derivatives in terms of tangents. This observation is useful when 
one wants to approximate an irregular function by a piecewise smooth (scaling) function. 

To illustrate the definitions of local fractional differentiability and critical order consider an 
example of a polynomial of degree n with its graph passing through the origin and for which the 
first derivative at the origin is not zero. Then all the local fractional derivatives of order less than 
or equal to one exist at the origin. Also all derivatives of integer order greater than one exist, as 
expected. But local derivatives of any other order, e.g. between 1 and 2 [see equations (||) and @] 
do not exist. Therefore critical order for this function at x = is one. In fact, except at a finite 
number of points where the function has a vanishing first derivative, critical order of a polynomial 
function will be one, since the linear term is expected to dominate near these points. 
Remark: We would like to point out that there is a multiplicity of definitions of a fractional deriva- 
tive. The use of a Riemann-Liouville definition, and other equivalent definitions such as Grunwald's 
definition, are suitable for our purpose. The other definitions of fractional derivatives which do not 
allow control over both the limits, such as Wyel's definition or definition using Fourier transforms, 
are not suitable since it would not be possible to retrieve the local nature of the differentiability 
property which is essential for the study of local behavior. Also, the important difference between 
our work and the work of 0, ^ is that while we are trying to study the local scaling behavior these 
works apply to asymptotic scaling properties. 



3 Fractional Differentiability of Weierstrass Function 

Consider a form of the Weierstrass function as given above, viz., 



W\(t) = ^ (s ~ 2)k sin \ h t, A > 1. (14) 
k=i 



Note that W x (0) = 0. Now 



d q W x (t) - (s _ 2)fc ^sin(A fc t) 

dtfl ' ^ dti 



provided the right hand side converges uniformly. Using, for < q < 1, 

d q sin(a;) d q ~ x cos(a;) 
dx q ~ dxi- 1 ' 
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we get 



dti d(XH) 



From the second mean value theorem it follows that the fractional integral of cos(A fc t) of order q — 1 
is bounded uniformly for all values of X k t. This implies that the series on the right hand side will 
converge uniformly for q < 2 — s, justifying our action of taking the fractional derivative operator 
inside the sum. 

Also as t — > for any k the fractional integral in the summation of equation fllED goes to zero. 
Therefore it is easy to see from this that 

E) q Wx(0) = lim dQWx ® =0 for q < 2 - s. (16) 

This shows that the gth local derivative of the Weierstrass function exists and is continuous, at 
t = 0, for q < 2 - s. 

To check the fractional differentiability at any other point, say r, we use t' = t — r and W(t') = 
W{t' + r) - W(t) so that W(0) = 0. We have 

oo oo 

W x (t') = ^A^^sinA^t' + ^-^A^-^sinAV 

k=l k=l 

oo 

= J2 x(s ~ 2)k ( cosXkrsinXkt ' + sinA fe r(cosAV - 1)). (17) 
fc=i 

Taking the fractional derivative of this with respect to t' and following the same procedure we can 
show that the fractional derivative of the Weierstrass function of order q < 2 — s exists at all points. 

For q > 2 — s, right hand side of the equation ( |I5|) seems to diverge. We now prove that the 
LFD of order q > 2 — s in fact does not exist. We do this by showing that there exists a sequence 
of points approaching along which the limit of the fractional derivative of order 2 — s < q < 1 
does not exist. We use the property of the Weierstrass function [|Hj], viz., for each t! G [0, 1] and 
< 5 < 5q there exists t such that \t — t'\ < 5 and 

c5 a < \W[t) - W(t')\, (18) 

where c > and a = 2 — s, provided A is sufficiently large. We consider the case of t' — and 
t > 0. Define g(t) = W{t) - ct a . 

Now the above mentioned property, along with continuity of the Weierstrass function assures 
us a sequence of points t± > t 2 > ... > t n > ... > such that t n — > as n — ► oo and g(t n ) = and 
g(t) > on (t n , e) for some e > 0, for all n (it is not ruled out that t n may be zero for finite n). 
Define 

g n (t) = 0, if t < t n , 

= g(t), otherwise. 

Now we have, for < a < q < 1, 

d"g n (t) 1 d g(y) 



d{t-t n y r(l -q)dtJ tn (t-y)i 



;dy, 
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where t n < t < t n -\. We assume that the left hand side of the above equation exists for if it does 
not then we have nothing to prove. Let 



tn (t ~ VY 

Now h(t n ) = and h(t n + e) > 0, for a suitable e, as the integrand is positive. Due to continuity 
there must exist an e' > and e' < e such that h(t) is increasing on (t n , e'). Therefore 

< f 9n{t \ \t-t , 71 = 1,2,3,.... (19) 
~ d(t-t n )i 1 v ; 



This implies that 



c rf9t " I <J!^L, „_12 3 f20) 



But we know from Eq. (^) that, when < a < q < 1, the left hand side in the above inequality 
approaches infinity as t — > 0. This implies that the right hand side of the above inequality does 
not exist as t — > 0. This argument can be generalized for all non-zero i' by changing the variable 
t" = t — t' . This concludes the proof. 

Therefore the critical order of the Weierstrass function will be 2 — s at all points. 



Remark: Schlesinger et al [25] have considered a Levy flight on a one dimensional periodic lattice 



where a particle jumps from one lattice site to other with the probability given by 

u — 1 00 

P(x) = ^— E u'Wx, +V) + S(x, -V)], (21) 

where x is magnitude of the jump, b is a lattice spacing and b > u > 1. 5(x, y) is a Kronecker delta. 
The characteristic function for P(x) is given by 

< i — 1 00 

P(fc) = E cos(Vk). (22) 

which is nothing but the Weierstrass cosine function. For this distribution the Levy index is 
logcu/logfe, which can be identified as the critical order of P(k). 

More generally for the Levy distribution with index \x the characteristic function is given by 

P(k) = Aexpc\k\^. (23) 

The critical order of this function at k = also turns out to be same as \i. Thus the Levy index 
can be identified as the critical order of the characteristic function at k = 0. 

4 Connection between critical order and the box dimension 
of the curve 

Theorem 1 Let f : [0, 1] — > M be a continuous function. 
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a) If 



*-*v [d(x — y)\ q 



where q,a G (0, 1), for all y G (0, 1), then dirriBf(x) < 2 — a. 
b) If there exists a sequence x n —> y as n — > oo snc/i i/iai 

gOM - f(y)) , f ^ 

llm U? \Ta = ±0 °' fOT 1 > ' 

[d(x„ - y)} q 

for all y, then dim^f > 2 — a. 

Proof: (a) Without loss of generality assume y = and /(0) = 0. We consider the case of q < a. 
As < q < 1 and /(0) = we can write J7| 



r( g )y (x-y)-^ 1 

Now 

As, by assumption, for q < a, 

lim^=0, 

a=->o ax 9 

we have, for any e > 0, a 5 > such that \d q f(x)/dx q \ < e for all x < 5, 

- T{q)l (x-y)-^ 
e 



dy. (24) 



-x q . 



r(g + i) 

As a result we have 

< for x<<5. 

Now this argument can be extended for general y simply by considering x — y instead of x and 
f(x) — f(y) instead of f(x). So finally we get for q < a 

\f{x)-f(y)\ < K\x-y\*, ioT\x-y\<6, (25) 

for all y G (0, 1). Hence we have [|i~0f 

dims/(x) < 2 — a. 
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b) Now we consider the case q > a. If we have 

lim = oo, (26) 

then for given M\ > and 5 > we can find positive integer N such that \x n \ < 5 and d q f(x n ) / dx q n > 
Mi for all n > N. Therefore by Eq. <M) 



> * r- * 



T(q)Jo (ar n -y)- ,+1 

r( g + i) » 

If we choose 5 = then we can say that there exists x < 5 such that 

f(x) > ki5 q . (27) 

If we have 

Um d*f(x n ) 

X n >0 dXn 

then for given M 2 > we can find a positive integer N such that d q f(x n ) / dx q n < —M 2 for all n > N. 
Therefore 



t( v . -M 2 /■*» cfy 



r(g) Jo (x n - y)-9+i 
-M 2 



™3 



r(g + i) 

Again if we write 5 = x^, there exists x < 5 such that 

f{x) < -k 2 6 q . (28) 
Therefore by ( p7|) and (|28|) there exists x < 5 such that, for g > a, 

\f(x)\ > K5 q . 

Again for any y G (0, 1) there exists x such that for q > a and \x — y\ < 5 

\f(x)-f(y)\ > k6 q . 



Hence we have 10 



dinis/(x) > 2 — a. 

Notice that part (a) of the theorem above is the generalization of the statement that C 1 functions 
are locally Lipschitz (hence their graphs have dimension 1) to the case when the function has a 
Holder type upper bound (hence their dimension is greater than one). 
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Here the function is required to have the same critical order throughout the interval. We can 
weaken this condition slightly. Since we are dealing with a box dimension which is finitely stable 
|I(J , we can allow a finite number of points having different critical order so that we can divide the 
set in finite parts having the same critical order in each part. 

The example of a polynomial of degree n having critical order one and dimension one is consistent 
with the above result, as we can divide the graph of the polynomial in a finite number of parts 
such that at each point in every part the critical order is one. Using the finite stability of the box 
dimension, the dimension of the whole curve will be one. 

We can also prove a partial converse of the above theorem. 

Theorem 2 Let f : [0, 1] — > M be a continuous function. 

a) Suppose 

\f(x)-f(y)\<c\x-y\ a , 
where c> 0, < a < 1 and \x — y\ < 5 for some 5 > 0. Then 

v dKfjx) - f(y)) 

hm — — = 0, for q < a, 

v^y [d{x-y)\i 

for all y G (0, 1). 

b) Suppose that for each y G (0, 1) and for each 5 > there exists x such that \x — y\ < 5 and 

\f(x)-f(y)\>c5 a J 

where c > 0, 5 < 5 for some 5 > and < a < 1. Then there exists a sequence x n — > y as 
n — » oo such that 

r d?(f(x n ) - f(y)) 

hm — -— = ±oo, for q > a, 

n-oo [d(x n -y))i 

for all y. 
Proof 

a) Assume that there exists a sequence x n — > y as n — > oo such that 

r d*(f(x n ) - f(y)) 

hm — t— = ±oo, tor q < a , 

™ [d{x n -y)]i 

for some y. Then by arguments between Eq. (^6|) and Eq. (p7[) of the second part of the previous 
theorem it is a contradiction. Therefore 



Now if 



hm — — rr = const or U, tor q < a. 

[d(x - y)]i 



lim — tCt—^ tt~~ = const, for q < a, 

x^y [d(x - y)Y 
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then we can write 

d q (f(x) - f(y)) 
[d{x - y)Y 



K + r){x,y), 



where K = const and rj(x, y) —>■ sufficiently fast as x — > y. Now taking the e derivative of both 
sides, for sufficiently small e we get 

d? + Hf(x) - f(y)) K(x-y)- e d^ix.y) r 
[d(x - r(l-e) [d(x-3/)]« 

As x — > y the right hand side of the above equation goes to infinity (the term involving rj does not 
matter since rj goes to sufficiently fast) which again is a contradiction. Hence the proof. 

b)The proof follows by the method used in the previous section to show that the fractional 
derivative of order greater than 2 — a of the Weierstrass function does not exist. 

These two theorems give an equivalence between the Holder exponent and the critical order of 
fractional differentiability. 

5 Local Fractional Derivative as a tool to study pointwise 
regularity of functions 

Motivation for studying pointwise behavior of irregular functions and its relevance in physical 
processes was given in the Introduction. There are several approaches to studying the pointwise 



behavior of functions. Recently wavelet transforms |29| , |30|j were used for this purpose and have met 
with some success. In this section we argue that LFDs is a tool that can be used to characterize 
irregular functions and has certain advantages over its counterpart using wavelet transforms in 
aspects explained below. Various authors |]32| , [40(1 have used the following general definition of 
Holder exponent. The Holder exponent ac(y) of a function / at y is defined as the largest exponent 
such that there exists a polynomial P n (x) of order n that satisfies 

\f(x)-P n (z-y)\=0{\x-y\«), (29) 

for x in the neighborhood of y. This definition is equivalent to equation (f25|), for < a < 1, the 
range of interest in this work. 

It is clear from theorem I that LFDs provide an algorithm to calculate Holder exponents and 
dimensions. It may be noted that since there is a clear change in behavior when order q of the 
derivative crosses the critical order of the function it should be easy to determine the Holder 



exponent numerically. Previous methods using autocorrelations for fractal signals |L0| involve an 
additional step of finding an autocorrelation. 



5.1 Isolated singularities and masked singularities 

Let us first consider the case of isolated singularities. We choose the simplest example f(x) = 
ax a , < a < 1, x > 0. The critical order at x = gives the order of singularity at that point 
whereas the value of the LFD iD 9=Q /(0), viz aT(a + 1), gives the strength of the singularity. 
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Using LFD we can detect a weaker singularity masked by a stronger singularity. As demonstrated 
below, we can estimate and subtract the contribution due to the stronger singularity from the 
function and find out the critical order of the remaining function. Consider, for example, the 
function 

f(x) = ax a + bx 13 , 0<a<P <1, x > 0. (30) 

The LFD of this function at x = of the order a is E) a f(0) = aT(a + 1). Using this estimate of 
stronger singularity we now write 

E> a f(0) 

G(x ] a) = f(x)-f(0)-—^-x a , 
which for the function / in Eq. (|30| ) is 

-x^ q . (31) 



d q G(x'a) bT((3+l) 



dxi T((3-q + l)' 

Therefore the critical order of the function G, at x = 0, is (3. Notice that the estimation of the 
weaker singularity was possible in the above calculation just because the LFD gave the coefficient of 
x a /T(a + 1). This suggests that using LFD, one should be able to extract the secondary singularity 
spectrum masked by the primary singularity spectrum of strong singularities. Hence one can gain 
more insight into the processes giving rise to irregular behavior. Also, one may note that this 
procedure can be used to detect singularities masked by regular polynomial behavior. In this way 
one can extend the present analysis beyond the range < a < 1, where a is a Holder exponent. 

A comparison of the two methods of studying pointwise behavior of functions, one using wavelets 
and the other using LFD, shows that characterization of Holder classes of functions using LFD is 
direct and involves fewer assumptions. The characterization of a Holder class of functions with 
oscillating singularity, e.g. f(x) = x a sin(l/x /3 ) (x > 0, < a < 1 and (3 > 0), using wavelets needs 
two exponents |3l|. Using LFD, owing to theorem I and II critical order directly gives the Holder 
exponent for such a function. 

It has been shown in the context of wavelet transforms that one can detect singularities masked 
by regular polynomial behavior |32j by choosing the analyzing wavelet with its first n (for suitable 
n) moments vanishing. If one has to extend the wavelet method for the unmasking of weaker 
singularities, one would then require analyzing wavelets with fractional moments vanishing. Notice 
that one may require this condition along with the condition on the first n moments. Further 
the class of functions to be analyzed is in general restricted in these analyses. These restrictions 
essentially arise from the asymptotic properties of the wavelets used. On the other hand, with the 
truly local nature of LFD one does not have to bother about the behavior of functions outside our 
range of interest. 

5.2 Treatment of multifractal function 



Multifractal measures have been the object of many investigations |34|, |36], p7[| . This formal- 



ism has met with many applications. Its importance also stems from the fact such measures are 



natural measures to be used in the analysis of many phenomenon []3~8 , j39fl . It may however happen 



that the object one wants to understand is a function (e.g., a fractal or multifractal signal) rather 
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than a set or a measure. For instance one would like to characterize the velocity of fully developed 
turbulence. We now proceed with the analysis of such multifractal functions using LFD. 

Now we consider the case of multifractal functions. Since LFD gives the local and pointwise 
behavior of the function, conclusions of theorem I will carry over even in the case of multifractal 
functions where we have different Holder exponents at different points. Multifractal functions have 
been defined by Jaffard [40] and Benzi et al. [IT|. However as noted by Benzi et al. their functions 



are random in nature and the pointwise behavior can not be studied. Since we are dealing with non- 
random functions in this paper, we shall consider a specific (but non-trivial) example of a function 
constructed by Jaffard to illustrate the procedure. This function is a solution F of the functional 
equation 

F(x)=J2\ l F(Sr\ x ))+g(x), (32) 
i=i 

where Si's are the afline transformations of the kind Si(x) = \i%x + 6, (with < 1 and fej's real) 
and Aj's are some real numbers and g is any sufficiently smooth function (g and its derivatives 
should have a fast decay). For the sake of illustration we choose H\ = /12 = 1/3, 61 = 0, b 2 = 2/3, 
Ai = 3- a , A 2 = 3-" (0 < a < p < 1) and 

g(x) = sin(27rx), if 16 [0,1], 
= 0, otherwise. 

Such functions are studied in detail in Ref. [f|(J using wavelet transforms where it has been shown 



that the above functional equation (with the parameters we have chosen) has a unique solution F 
and at any point F either has Holder exponents ranging from a to (3 or is smooth. A sequence of 
points 5^(0), Si 2 S'i 1 (0), • • ■ , S{ n ■ ■ ■ (0) , ■ ■ •, where ik takes values 1 or 2, tends to a point 
in [0, 1] (in fact to a point of a triadic cantor set) and for the values of /i^s we have chosen this 
correspondence between sequences and limits is one to one. The solution of the above functional 
equation is given by Ref. |4(| as 



i^) = E E K---Kg(s- 1 ---sr 1 1 (x)). (33) 

n=0 »!,•••) in=l 



Note that with the above choice of parameters the inner sum in (|33|) reduces to a single term. 
Jaffard has shown that 

h(y) = hminf |° g £*M \ " " ^ (34) 

where {h(y) ■ ■ ■ i n (y)} is a sequence of integers appearing in the sum in equation (^) at a point y, 
and is the local Holder exponent at y. It is clear that h m i n = a and h max = j3. The function F at 
the points of a triadic cantor set have h(x) G [at, /3] and at other points it is smooth ( where F is as 
smooth as g). Now 

d«(F(x) - F(y)) 



2 

l d ( x - y)\ q ~ S n,5=i Ail '" Al " 
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di[g(S- 1 ---S-\x))-g(S^...Sr\y))} 
[d(x - y)Y 

oo 2 

= 12 12 K ■ ■ ■ K(vh ■ ■ ■ Vi„)~ q 

n=0 h,—,i n =l 

dojgjS- 1 ■ ■ ■ Sr\x)) - g(S£ ■ ■ ■ Sr\y))] 

MS' 1 ■ ■ ■ S~\x - y))]« ' [6b) 

provided the RHS is uniformly bounded. Following the procedure described in section III the 
fractional derivative on the RHS can easily be seen to be uniformly bounded and the series is 
convergent if q < mm{h(x), h(y)}. Further it vanishes in the limit as x — > y. Therefore if q < h(y) 
E) q F(y) = 0, as in the case of the Weierstrass function, showing that h(y) is a lower bound on the 
critical order. 

The procedure of finding an upper bound is technical and lengthy. It is carried out in the 
Appendix below. 

In this way an intricate analysis of finding out the lower bound on the Holder exponent has been 
replaced by a calculation involving few steps. This calculation can easily be generalized for more 
general functions g{x). Summarizing, the LFD enables one to calculate the local Holder exponent 
even for the case of multifractal functions. This fact, proved in theorems I and II is demonstrated 
with a concrete illustration. 



6 Conclusion 

In this paper we have introduced the notion of a local fractional derivative using Riemann-Liouville 
formulation (or equivalents such as Grunwald's) of fractional calculus. This definition was found to 
appear naturally in the Taylor expansion (with a remainder) of functions and thus is suitable for 
approximating scaling functions. In particular we have pointed out a possibility of replacing the 
notion of a tangent as an equivalence class of curves passing through the same point and having the 
same derivative with a more general one. This more general notion is in terms of an equivalence 
class of curves passing through the same point and having the same critical order and the same 
LFD. This generalization has the advantage of being applicable to non-differentiable functions also. 

We have established that (for sufficiently large A ) the critical order of the Weierstrass function 
is related to the box dimension of its graph. If the dimension of the graph of such a function is 
1 + 7, the critical order is 1 — 7. When 7 approaches unity the function becomes more and more 
irregular and local fractional differentiability is lost accordingly. Thus there is a direct quantitative 
connection between the dimension of the graph and the fractional differentiability property of the 
function. This is one of the main conclusions of the present work. A consequence of our result is 
that a classification of continuous paths (e.g., fractional Brownian paths) or functions according to 
local fractional differentiability properties is also a classification according to dimensions (or Holder 
exponents). 

Also the Levy index of a Levy flight on a one dimensional lattice is identified as the critical order 
of the characteristic function of the walk. More generally, the Levy index of a Levy distribution is 
identified as the critical order of its characteristic function at the origin. 

We have argued and demonstrated that LFDs are useful for studying isolated singularities and 
singularities masked by the stronger singularity (not just by regular behavior). We have further 
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shown that the pointwise behavior of irregular, fractal or multifractal functions can be studied using 
the methods of this paper. 

We hope that future study in this direction will make random irregular functions as well as 
multivariable irregular functions amenable to analytic treatment, which is badly needed at this 
juncture. Work is in progress in this direction. 
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Appendix: Upper bound on critical order of function of 
Eq. gOj 



Our aim in this is appendix is to find an upper bound on the critical order and hence on the local 
Holder exponent for the function of equation (p3|). Our procedure will be similar to that of Jaffard 
40fl the only difference being that we take the fractional derivative instead of the wavelet transform. 



We proceed as follows. The defining equation for F(x) when iterated N times gives 

N-l 2 

H*) = E E K---\ in g(sr n l ---sr l 1 (x)) 

71=0 H,---,j n = l 

+ E A n ---A, JV F(Sr;...5r 1 (x)), are [0,1]. (36) 

We now consider 

df(F(x) - F(y)) ™ ^ , w __ 



[d(x - y)]i 



E E ■ " ^ini^h " '/O 

n=0 ii, •••,j 7l =l 

d^g(S^...Sr\ x ))-g(SrJ...Sr l (y))} 
[d(Sri ■ ■ ■ Sr\x - y))]i 

2 

+ E " ' ^*Jv(A*m. " ' Mijv) 9 

ijv=l 

g[FQSg ■ ■ ■ Sr\x)) - F(S£ ■ ■ ■ Sr\y))} 
[d{S£ - - Sr\ x - y))]* 



h, — ,iw=l 



(37) 



Let us denote the first term on the RHS by A and the second term by B. In the following y e (0, 1). 
Choose N such that 3^( iV + 1 ) < \x — y\ < 3~ N . Denote \ n (y) — Ki(y) ' ' '\ n {y)- F° r the values of /XjS 
we have chosen fi^ • ■ ■ = 3~ n . Now since g is smooth \g(x) — g(y)\ < C\x — y\, 

d<1 [9(S~! ■ ■ ■ Sr l ( x )) - gjSr 1 ■ ■ ■ S-\y))] x-y x _ q 
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From the way we chose /i^s and \x — y\ this term is bounded by C3 n ^ N ^ q \ Therefore the 
first term in equation (|37D above is bounded by 



N-l 



A < \n(y)T g 3 n{1 - g) 3- N{1 - q) 

n=0 

N-l 

= OS-* 1 -* £ 

v - x n 4- 

3A(Ar-i)(y) 3 2 A(at_i)(^) 



n=0 

AfAr-ov,/! Ac 



< C3- JV ( 1 ^A (jV _ 1)fa) 3 JV - 1 (l + o y- 2)fa) + o y- 3)(y) +■■■)• (39) 



The quantity in the brackets is bounded. Therefore 

A < C3 N «\ N{y) . (40) 



Now we consider the second term on the RHS of equation (|37|) and find a lower bound on that 
term. We assume that F(x) is not Lipschitz at y (for otherwise the LHS of equation ( 1371) will be 
zero in the limit x — > y and the case is uninteresting). Therefore there exists a sequence of points 
x n approaching y such that 

\F{x n ) - F{y)\ > c\x n - y\. (41) 

That the function is not Lipschitz at y implies that it not Lipschitz at ■ ■ ■ (y). Therefore we 
get (one may recall that only one of the 2 n terms in the summation ... i N =\ is non-zero) 

B > c\ N{y) 3 N<1 3 N ^3- N ^ 

> c\ N{y) 3 N(1 . (42) 

Therefore there exists a sequence x n such that 

Since (p[) is valid for every c for large enough n, 

d q (F(x n ) — F(y)) , Nq 
[d{x n -y)Y ^ CXn ^ ' 

where C = c — C . Now from the expression of h(y) it is clear that when q > h(y) LFD does not 
exist and the critical order is bounded from above by h(y). 
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